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A consequence of strong antiferromagnetic correlations in 
models of high-Tc cuprates is the appearance in photoemis- 
sion (PES) calculations of considerable more weight above the 
Fermi momentum pp than expected for non-interacting elec- 
trons. This effect, first discussed by Kampf and Schrieffer 
(Phys. Rev. B 41, 6399 (1990)) under the name of "shadow 
bands", is here analyzed in the two dimensional Hubbard 
and t-J models using Monte Carlo and exact diagonalization 
techniques. It is concluded that for underdoped compounds 
weight above pp could be observable in experimental PES 
data, while in the overdoped regime it will be likely hidden in 
the experimental background. In the intermediate region the 
signal is weak. Our results are thus compatible with recent 
experimental data by Aebi et al. (Phys. Rev. Lett. 72, 2757 
(1994)) for Bi2212 at optimal doping. However, to definitely 
prove the existence of this effect in the cuprates, experiments 
in the underdoped regime are necessary. 



The importance of antiferromagnetic correlations in 
the normal and superconducting states of the high criti- 
cal temperature cuprate materials is under much discus- 
sion. While recently most of the debate has been con- 
centrated on the symmetry of the superconducting order 
parameter, studies of the strength of the antiferromag- 
netic correlation length, ^af , in the normal state are still 
crucial to test these ideas. A key issue is how large £,af 
should be to produce observable effects in experiments 
for the high-Tc compounds. NMR studies in the normal 
state of optimally doped YBa2Cu306.94 (YBCO) suggest 
S,AF ^ 2a (where a is the lattice spacing). H] Naively, 
this correlation seems too small to be of relevance. On 
the other hand, recent photoemission (PES) experimen- 
tal results by Aebi et al. @ in BiaSraCaCuaOg (Bi2212) 
with Tc=85K, using sequential angle-scanning data ac- 
quisition to obtain PES intensities within a narrow en- 
ergy window near the Fermi energy Ep, reported evi- 
dence compatible with the antiferromagnetically induced 
"shadow bands" first predicted by Kampf and Schrieffer. 
PI At half-filling, these bands are caused by the enlarged 
magnetic unit cell of the Cu02 planes produced by the 
long range antiferromagnetic order in the ground state. 
This effective reduction in the size of the Brillouin zone 
(_BZ) has interesting implications for PES experiments. 
Q For example, along the diagonal px = Py = p, and as- 
suming long-range order, peaks at momenta pi = (p, p) 
and p2 = (tt — p, TT — p) should appear at the same en- 
ergy location, for any value of p. The PES weight ob- 



served in the region above the non-interacting pp is in- 
duced by strong magnetic correlations. |3|] 

How does this antiferromagnetically induced PES sig- 
nal evolve as ^af is made finite by the effect of hole 
doping? It is likely that its intensity will smoothly di- 
minish doping away from half-filling. Since in optimally 
doped Bi2212 the antiferromagnetic correlation length 
in the Cu02 planes should be similar to that of opti- 
mally doped YBCO (since they should have the same 
in-plane hole density), then a natural question arises: 
can a small correlation length (~ 2a) produce observable 
weight in PES experiments above pp? On one hand, re- 
cent calculations ^ carried out in an antiferromagnetic 
background have reproduced the flat band features near 
p = (tt, 0) observed in Bi2212 suggesting that a short S^af 
can appreciably influence some experimental quantities. 
On the other hand, since the actual experimental PES 
signal (Fig. la) is weak, concerns may arise about the 
interpretation of the data (as a reference, PES weight 
induced by antiferromagnetism in recent experiments for 
the insulating layered copper oxide Sr2Cu02Cl2 is shown 
in Fig. lb. [g| ). Thus, a theoretical calculation is needed 
to compare with experiments, and decide if Aebi et al.'s 
PES data are compatible with models of strongly corre- 
lated electrons which present PES weight above pp. 

Here, this issue is explicitly addressed. PES spectra 
and spin correlations are calculated for models of cor- 
related electrons expected to qualitatively describe the 
Cu02 planes. Consider flrst the two dimensional (2D) 
one band Hubbard model defined by the Hamiltonian 
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in the standard notation. We simulated this model nu- 
merically using standard quantum Monte Carlo tech- 
niques. In order to extract the dynamical spectral func- 
tion A(p,a;) corresponding to the removal or addition of 
an electron with momentum p to the system, the max- 
imum entropy (ME) technique was used. P| To analyze 
the strength of the signal above pp, we calculated the 
amount of spectral weight below the chemical potential 
/i at momenta along the diagonal px — Py in the BZ with 
respect to the total intensity (adding PES and inverse 
PES) which for the Hubbard model satisfies the sum rule 
/_ dujA{p, w) = 1 at all dopings. 

Before describing the computational results, let us clar- 
ify when a theoretically calculated PES signal can be la- 
beled as being "observable" in an experiment. PES spec- 
tra have large backgrounds, whose origin and shape are a 



matter of discussion, superimposed on the actual relevant 
signal. This background depends on the momentum, and 
also changes from sample to sample with fluctuations as 
large as 50%. Since the background is convex, the natural 
requirement for a PES theoretical signal to be observable 
is that the combination background-signal produces a lo- 
cal maximum (i.e. a peak in the measured intensity). R] 
From the data shown in Fig. lb, and the intensity of the 
signal at the last point where the dispersion is observed 
i.e. p = (O.Ttt, O.Ttt), we believe that a peak with an 
intensity of roughly about 10% of the largest signal (lo- 
cated at p = (O.Stt, O.Stt)) would be at the verge of being 
detected. R] This is the criterion followed in the present 
paper to label a result as "observable" . 

The choice of coupling is important in our search for 
PES weight above pp. For example, we observed that 
working on an 8 x 8 cluster, at U/t = 4, half-filling and 
temperature T = t/4, the percentage of PES spectral 
weight at p = (37r/4, 37r/4), i.e. the next available mo- 
mentum after (7r/2, 7r/2) on this cluster, is very small 
(less than 5% of the total), even though the spin corre- 
lations show clear indications of long-range order. Then, 
the actual value of the local moments is as important as 
the antiferromagnetic correlation length for the effect we 
are investigating. Since evidence has been recently given 
that another feature induced by antiferromagnetism, |^] 
namely the "hole" pockets, are washed out by temper- 
ature effects in QMC simulations at U/t = 4, then we 
consider this coupling to be too small for our purposes. 
Actually, studies of the optical conductivity have shown 
that a larger coupling U/t, approximately between 8 and 
12, is needed to correctly reproduce the main features 
observed experimentally, [pj Since at U/t = 12 there are 
serious numerical instabilities in the simulations, then 
here the analysis was restricted to U/t = 8. 

In Fig.2, A(p, w) at T = t/2 is shown. ^ It is difficuh 
to reduce T due to sign problems, but nevertheless this 
temperature allows us to study the PES signal above pp 
at different correlation lengths as the density is changed, 
which is the main purpose of the paper. At half-filling, 
(n) — 1, the chemical potential is located in the gap. The 
percentage of spectral weight is shown for each momen- 
tum. A nonzero signal above the non-interacting Fermi 
momentum is clearly visible, and at p = (37r/4, 37r/4) it 
carries ~ 23% of the total weight. This result is very 
similar if the temperature is reduced to T = t/4, thus fi- 
nite temperature effects are not too severe for this quan- 
tity at half-filling. Actually, our results for the inten- 
sity of the weight above pp are in excellent agreement 
with the spin density wave mean-field approximation. |ll| ] 
Since for a pure spin- 1/2 antiferromagnet the weight at 
p = (7r/4, 7r/4) and p = (37r/4, 37r/4) should be identical 
at this density, it is natural to conclude that the finite 
coupling U/t is responsible for the reduction of the inten- 
sity of the PES signal above the naive pp at half-filling. 

Away from half-filling, at (n) = 0.87, the amount of 
weight at p = (37r/4, 37r/4) is reduced to ~ 10%, which 
is still visible in the scale of the plot. The height 



of the peak, as a percentage of the peak height at 
p = (7r/2, 7r/2) and half-filling is about 15%. Follow- 
ing the criteria described before, the result obtained at 
(n) — 0.87 in the Hubbard model is labeled as still "ob- 
servable" (although it is rather weak). Whether this 
PES weight corresponds to an actual sharp band dis- 
persing, to a broader feature, or a combination of both, 
is difficult to address with the ME technique which has 
low resolution, but nevertheless it is clear that it is in- 
duced by antiferromagnetism. At (n) = 0.70, the sig- 
nal at p = (37r/4, 37r/4) carries a small weight of only 
4%, and as the electronic density is reduced further the 
system smoothly converges to the non- interacting limit. 
Then, for this particular calculation we tentatively con- 
clude that a doping of 25% holes makes the antiferromag- 
netically induced weight almost negligible, while at 12% 
doping the effect is still observable. To relate these results 
with experiments, in Fig. 3a we show the spin correlations 
calculated numerically. At (n) = 1, the correlation is 
robust but it decays slowly to zero due to temperature 
effects, while at (n) = 0.70 it is clearly very small. At an 
intermediate density (n) — 0.87, ^af is between one and 
two lattice spacings, resembling the experimental situa- 
tion in YBCO, and presumably also in Bi2212, since both 
are at optimal doping. ||l2] In spite of this fact, the weight 
above pp is observable in Fig.2, showing that the effect 
is mainly caused by the local short-distance antiferro- 
magnetic spin arrangement. Thus, the presence of large 
spin correlation lengths is not a necessary condition for 
the observation of weight above the noninteracting Fermi 
momentum. 

To check the sensitivity of our conclusions to the model 
used, let us also consider the well-known 2D t — J Hamil- 
tonian. This model cannot be studied with Monte Carlo 
techniques, thus exact diagonalization was used. M This 
algorithm works at T=0, and dynamical information 
can be obtained in real time. Its restriction to small 
clusters should not be a major problem in calculations 
where ^af is very small. To increase the momentum 
resolution along the diagonal in the BZ, we combined 
the results of the 16 sites cluster (providing momenta 
(0, 0), {tt/2, 7r/2), (tt, tt)) and the 18 sites cluster (contain- 
ing(V3,7r/3),(2V3,27r/3)). |1 

In Fig. 4, the PES A(p, ui) spectrum is shown for the 
t-J model. J/t = 0.4 was selected as an example, but 
we checked that the results are qualitatively similar in 
the range between J/t = 0.2 and J/t = 0.8. As expected, 
at half-filling the largest peak near the chemical poten- 
tial (quasiparticle) is obtained at p = (7r/2,7r/2). In- 
creasing the diagonal momenta away from it, a con- 
siderable amount of spectral weight induced by S,af is 
observed. Moving away from half-filling into the sub- 
space of two holes (nominal density (n) ^ 0.88) the 
dominant peak remains at p = (7r/2,7r/2) within our 
momentum resolution. At (7r/3, 7r/3) the quasiparticle 
strength is still large and coherent. On the other hand, 
at p = (27r/3, 27r/3) the peak seems now broader in the 
scale used, although its integrated spectral weight re- 



mains close to that of p = (tt/S, tt/S). The height of the 
peak at p = (27r/3, 27r/3) as a percentage of the largest 
peak located at p = (7r/2,7r/2) (with or without holes) 
is 15-20% i.e. within the "observable" region defined be- 
fore. Finally, at density (n) = 0.77, the result resembles 
that of a non-interacting system with a Fermi momen- 
tum close to p = (7r/3,7r/3), above which the signal is 
too weak to be observable in PES experiments. Then, 
our rough estimations within the t-J model are similar 
to those of the Hubbard model, i.e. weight above pp 
can still be observed at (n) ~ 0.88 but no longer at den- 
sity (n) ~ 0.77. To make contact with experiments it is 
again necessary to consider the spin correlations which 
are shown in Fig. 3b. At half-filling, ^af is clearly larger 
than the lattice size. At (n) ~ 0.88, a crude exponential 
fit of the spin correlation vs. distance gives ^af ~ l-5a 
(similar to that of YBCO and Bi2212 at optimal doping), 
while at (n) ~ 0.77, ^af is less than one lattice spacing. 
Then, we again arrive to the conclusion that for a real 
material with ^af ~ 2a the antiferromagnetically gen- 
erated PES weight is weak but may still be observable 
above the background. 

In Fig. 5, A(p,w) is shown again at (n) ~ 0.88 using 
an expanded energy scale. The dispersion of the sharp 
peak (I) discussed before in Fig. 4, has a bandwidth of 
order J, while at higher energies a considerable amount 
of spectral weight is found contributing to the bulk of the 
valence band (II). This double peak structure is in good 
agreement with the prediction of Kampf and Schrieffer 0] 
i.e. the lowest energy peak (I) and the incoherent part 
(II) at p = (27r/3, 27r/3) are the shadow bands caused by 
antiferromagnetic correlations. ME entropy techniques 
cannot resolve these two peak structure, and thus exact 
diagonalization is needed to quantitatively analyze these 
features. 

Summarizing, an analysis of the PES spectra in the 2D 
Hubbard and t-J models, at several densities and cou- 
plings, was reported. If these models qualitatively repro- 
duce the physics of the high-Tc compounds, then we con- 
clude that the effect predicted by Kampf and Schrieffer 
based on a robust ^af should be observable even for ma- 
terials with spin correlations lengths of only a couple of 
lattice spacings, as in Bi2212 at optimal doping. This is 
compatible with the experimental results of Fig. la. How- 
ever, this regime is at the verge of observability. The 
PES signal above pf should no longer be visible above 
the large PES experimental background at slightly larger 
dopings. To gather further evidence that the weak exper- 
imental signal is indeed caused by antiferromagnetism we 
believe that it is necessary to carry out PES experiments 
as a function of hole doping. The strength of the signal 
above pp should increase as the system moves away from 
the optimal doping towards half-filling. A possible candi- 
date for such a study is YBCO with a critical temperature 
of about 60K. Another alternative within the Bi2212 fam- 
ily would be to consider Bi2Sr2Cai_xLuxCu208+5 and 
Bi2Sr2-xLaxCaCu208+5 which seem to be underdoped. 
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Figure Captions 

1. (a) PES intensity in Bi2212, as reported by Oster- 
walder et al. (Ref. |H|). Their method produces 
PES intensity at constant energy for aU momenta, 
while conventional methods provide complete PES 
energy distribution curves at a few locations in the 
BZ. Each solid line corresponds to a fixed energy 
scan starting at the bottom at 0.3 eV above Ep, 
and arriving at the last top line at Ep. The spectra 
were measured at a polar angle of 39'^, and for az- 
imuthal angles spaced l" apart beginning near the 
FM hne and ending near the FX hne. The "5 x 1" 
band is explained in the original text Ref. |0| and it 
is of no concern to us. We thank P. Aebi and J. Os- 
terwalder for providing us with these unpublished 
data; (b) PES intensity as a function of energy for 
SraCuOaCla taken from B. O. Wells et al. (Ref. H). 



The momenta are given in units of tt and along the 
diagonal in the 2D square lattice convention. 

2. A(p, Lu), evaluated using QMC and ME techniques, 
for the 2D Hubbard model at U/t = 8, T = t/2 on 
an 8 X 8 cluster, at several densities (n) . The mo- 
mentum label varies along the diagonal in the BZ 
in units of 7r/4, and the percentages correspond to 
the integrated PES part of the spectral weight with 
respect to the total intensity (=1). The energy is 
in units of t. 

3. (a) Spin-spin correlation 4(S?S?_^j)(-l)'Jl vs dis- 
tance, j = Ijj, for the 2D Hubbard model calculated 
using QMC at T = t/2, U/t = 8, and several densi- 
ties on an 8 X 8 cluster; (b) Spin-spin correlations as 
defined in Fig. 3a, for the 2D t-J model calculated 
using exact diagonalization techniques on a 4 x 4 
cluster with two holes, at several couplings. 

4. PES A(p,w) evaluated using exact diagonalization 
techniques for the 2D t-J model, at J/t = 0.4 on 4x 
4 and -s/lS x \/l8 clusters. The densities are shown 
in the figure. We assumed t = 0.4eV, and provided 
a width d = O.lt to the peaks. The momenta are 
indicated, and the relevant peaks are shaded. 

5. PES A(p,w) for the t-J model at (n) - 0.88, 
J/t — 0.4, clusters of 16 and 18 sites, and expanding 
the energy scale to observe the two peak structure. 
We use S = 0.25t and t = 0.4eV. 
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